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Abstract 


This paper deals with the design of multi-axis force (also known as 
force/torque) sensors, as considered within the framework of optimal 
design theory. Optimal design procedures consist of finding the com¬ 
bination of design variables that extremizes some optimality criterion: 
provided a suitable mathematical formulation of the problem, solutions 
can be efficiently obtained through currently available numerical tech¬ 
niques. The principal goal of this paper is to identify a mathematical 
objective function, whose minimization corresponds to the optimiza¬ 
tion of sensor accuracy. The methodology employed is derived from 
linear algebra and analysis of numerical stability. An objective func¬ 
tion which can be applied to a large class of sensor configurations is 
proposed. The problem of optimizing the number of basic transducers 
employed in a multi-component sensor is also addressed. Finally, ap¬ 
plications of the proposed method to the design of a simple sensor as 
well as to the optimization of a novel, 6-axis miniaturized sensor are 
discussed. 
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1 Introduction 


The development of multi-axis force sensors, i.e. instruments for measuring several (up to 
6) components of force and torque simultaneously, has been initially undertaken in fields 
such as wind-tunnel testing, adaptive control of machines and thrust stand testing of rocket 
engines. Typical in those areas is the problem of monitoring forces of variable directions 
and intensity. Some of the most interesting sensors designed in this phase are reviewed by 
Doebelin [1]. Starting from the mid seventies, a major impulse to research on such sensors 
has been given by the recognition of their large relevance to robotic and telemanipulation 
applications. 

The ability to sense the arm-environment interactions is a crucial need for robots to evolve 
from purely repetitive behaviors to some degree of autonomy in unstructured surroundings 
or tasks. The feed-back of forces and torques exerted by the arm’s end effector is not only 
instrumental for the accomplishment of most tasks involving modification of the environment 
by contact, but is also critical in guaranteeing safe operation of the arm. 

Considering the most common case of a serial link manipulator, a complete characteriza¬ 
tion of the system of forces acting on a portion of the arm can be obtained by force sensors 
interposed between that portion and the rest of the arm. Force sensing at intermediate 
sections of the kinematic chain is likely to be very useful for next-generation robots using 
the whole arm surface to interact with the environment [2, 3]. However, since the measure¬ 
ment of interaction forces can be disturbed by inertial forces arising from acceleration of the 
masses between the sensor and the end-effector, force sensors are usually placed as close as 
possible to the distal end of the arm, most often at the wrist. Several configurations have 
been proposed for such sensorized wrists (e.g. [4, 5, 6]), and some are commercially available. 
A pedestal force sensor has also been described in [4]. 

The increasingly demanding tasks assigned to automatic manipulation require ever more 
stringent accuracy, stiffness, encumbrance, weight, speed of response, ease of interfacing and 
reliability of force sensors for robotics. A particularly promising application of force sensing 
to fine manipulation, for instance, is the so-called intrinsic contact sensing concept [7, 8]. 
Very detailed information about the contact between two surfaces can be obtained based on 
force measurements and geometric considerations. By integrating force sensors in the very 
“fingertips” of the arm, interactions with manipulated objects can be monitored very closely. 
The improvement in manipulation dexterity has proven considerable for both simple parallel 
jaw grippers [9], and for dextrous articulated hands [10, 11, 12]. A drastic reduction of size 
and weight of the force sensors with respect to existing designs turned out to be mandatory 
for this application. 

The stringent and conflicting requirements of this and other advanced applications of 
force sensing render unsatisfactory the traditional, intuitive approach to force sensor design, 
and motivated the investigation of a more systematic approach. In this paper the design 
of multi-axis force sensors is considered within the framework of optimal design theory, 
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the branch of operations research currently offering the most developed tool for systematic 
synthesis of optimized designs. 


2 Problem Formulation 

According to Vanderplaats [13], an optimal design problem (which is basically a nonlinear 
programming problem) can be stated as follows: 

Find the set of r design variables Xi that minimizes the objective function F(x), 
subject to: 

1. p inequality constraints: 9i{x) < 0, 1 < i < p; 

2. q absolute constraints: hj(x) = 0, 1 < j < q', 

3. n bounds defining the feasible design region: Xk,i n f < < %k,sup, 1 < k < n. 

The variables characterizing a force sensor design can be chosen within a large range of 
quantities related to its geometric description, dimensions, material properties, machining 
process, transduction principle, etc.. Any quantity contributing to the identification of a 
design, which is not allowed to vary in the optimization process, is considered a design 
parameter. 

Inequality conditions are unilateral constraints that must be satisfied in order for the 
design to be acceptable. For example, stresses in structure members must not exceed specified 
values. The maximum displacement under nominal load, and the resolution of the sensor 
(intended here as the the minimum detectable load) can be conveniently considered under 
the form of inequality constraints. 

By means of an equality constraint, precise conditions that must be met by the design 
can be stated. For instance, these constraints can be used to define the system interfaces. A 
complex problem can be reduced to a simpler sub-optimal one by using equality constraints 
to assign some variables a fixed value. 

Finally, boundary constraints defining the feasible region for design variables prevent 
the optimization algorithm from converging to unacceptable solutions (e.g., sensor structure 
members with negative thickness). 

The objective function F(x) is the measure adopted for quantifying the quality of a 
design. A most important decision is implied by the choice of the quality to be optimized. 
In some cases emphasis can be put on cost, encumbrance, or weight reduction. Very often, 
however, the accuracy of measurements is of paramount importance in designing a force 
sensor. It is possible to take into account these and other possible criteria at the same time 
by using multi-criteria optimization techniques [14], e.g. by assigning weighting factors to 
each criterion. Although greater flexibility of the optimization tool can be achieved in this 
way, there is a risk with multi-criteria techniques of loosing insight into the design process. In 
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this paper, all the above mentioned and other possible requirements are considered as design 
constraints, and the focus is on the choice of an objective function capable of extracting, 
from the class of acceptable designs, the one giving best results in terms of measurement 
accuracy. In order to obtain general solutions from the optimal design method, we propose 
an objective function that covers a wide scope of design variables. In particular, an attempt 
is made at addressing the choice of the number of basic transducers to be employed in a 
multi-component sensor. 


3 Mathematical Model of a Force Sensor 

In general, a multi-axis force sensor is a device in which several simple transducers measure 
the effects of unknown loads on a mechanical structure. These basic measurements are 
processed in order to evaluate the components of the applied load. As far as a linear behavior 
can be hypothesized for the sensor, a model can be written as 


v = Cp (1) 

where p G 9? n , n < 6, is the unknown load vector, v € 5i m , m > n is a vector collecting 
the basic measurements, and C G 9? mXn is a constant, full rank matrix characteristic of the 
sensor. The unknown vector p is composed of some or all the components of the resultant of 
the load; each component is normalized with respect to its maximum value, which is given 
as a design specification. When basic measurements are relative to strains or displacements 
(as it most often occurs in practical applications), the C matrix is usually referred to as the 
compliance matrix of the sensor. 

The behavior assumed in equation 1 is ideal under several regards. For instance, the 
linearity assumption is rather strong. However, it should be pointed out that we are forced 
to such approximation by the need for a viable algorithm for inverting the model, i.e. for 
solving equation 1. Nonlinear elastic models of even simple mechanical structures are too 
complex for being inverted, and even though a tabulated calibration of a sensor is conceivable, 
there is no practical application of such technique to force sensors. As a matter of fact, linear 
(or piece-wise linear) approximations can be sharpened by calibration techniques (see e.g. 

[15]) - 

As a second remark about the model 1, we note that the measurement vector v will be 
in all likelihood affected by an error Sv due to inaccuracies of basic transducers. Similarly, 
also the imprecise knowledge of the actual C matrix, which will be affected by a (matrix) 
error SC, should be taken into account. Thus, a model for a real sensor should be written as 

v + 5v = (C + £C)(p + Sp), (2) 

where £p is the error resulting on the ultimate information from the sensor, the load vector 

P- 
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Consider for instance the origin of errors Sv and SC for a strain gauge based sensor. 
As is well known, strain gauges are electrical resistors that can be bonded on a deformable 
structure, and vary their resistance according to the mechanical strain of the structure. 
Errors 8v in evaluation of structure strains can be caused by electrical noise and thermal 
drift in resistance measurements, and by discretization errors in conversion of data from 
analog to digital. The elements of the compliance matrix C are in turn determined only 
approximately, either using elasticity theory formulas or a direct calibration procedure. Both 
these methods result in a relative error e c (usually much larger in the first case), arising 
from inappropriate modeling of sensor structure (for calculated C’s), or from inaccuracies 
in calibration. The effects of nonlinearities in the stress-strain relationship of the structure, 
the imperfect stiffness of the glue used to bond the gauges, and many other uncontrollable 
factors will also add to e c . 

Multi-axis force sensors can be subdivided between those using the minimal number of 
basic transducers for the measure of the unknown load components (i.e., vn — n), and those 
having redundant sensors (m > 7i). In the latter case, no exact solution of equation 1 is 
possible, and diverse approximate solutions can be chosen. This classification will be used 
in the following discussion. 


4 Minimal Sensors 

For m = n, the generalized form of Wilkinson’s formula for error propagation in linear 
algebraic systems (see [16]) can be utilized in order to give an a priori estimate of the 
relative error on p: 


e p - (e v + e c ) K p ( C), (3) 

where = M, e c = e p = are respectively the relative errors on strain measure¬ 

ments, on calibration and on the results. The symbol || • || represents a generic vector norm, 
or the subordinate matrix norm if the argument is a matrix *. The propagation factor K p 
has the upper bound: 


K P { C) < 


N c 

1-iW 


where Nc is the condition number of the compliance matrix C, defined as 


(4) 


Nc = |iC|| lie- 1 !!, 


and where it is assumed Nc < 1. 

1 Vectors are designed by lower case boldface letters in this paper, while matrices are upper case boldface 
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A block scheme of equation 3 is given in figure 1. It appears from equation 4 that large 
condition numbers of C can spoil even the most accurately measured and calibrated sensor. 

After its definition, N c > 1, which implies K p > 1: this can be seen as an instance of the 
general principle of information theory regarding the entropy increase due to elaboration of 
data [17]. In the best case, when N c = 1 (which occurs only for orthogonal matrices and 
their multiplies), and K p = (1 — e c ) _1 ~ 1, the global error e p is simply the sum of source 
errors e v and e c . 

Given a sensor design, substantial reductions of e v and e c can be achieved only by using 
more sophisticated technologies, materials and components in the construction of the sensor, 
and finer models of the structure, or more accurate instrumentation for strain measurement 
and calibration. An improvement of these factors with respect to present sensors is therefore 
achievable by increasing their cost. Also assuming no cost constraint to the design, how¬ 
ever, there are absolute upper bounds to possible reduction of source errors set by present 
technological state-of-art, and by inherent measurement accuracy limitations. 

The amplification factor K p in equation 3 is related to the elements of the compliance 
matrix C, which can be chosen by the designer by varying, for instance, the thickness of 
some members in the structure, the position of strain-gauges, etc. Thus, the problem of 
optimizing the accuracy of multi-axis force sensors can be split in a technology-dominated 
sub-problem, and a design-dominated one. 

It should be pointed out that equation 3 does not assume any specific algorithm for 
the solution of the linear system 1. From this viewpoint, different algorithms such as LU 
factorization, singular value decomposition (SVD), and even special solution schemes relying 
on a possibly structured C matrix, are equivalent. The application of these methods with 
finite precision machines will produce numerical errors that are typically some orders of 
magnitude smaller than source errors, and can be considered as a small portion of the 
overall calibration error e c [18]. 

Since the right-hand member of equation 4 is a monotonic non-decreasing function of Nc, 
it seems natural to choose the condition number of the compliance matrix C as the objective 
function of the optimization procedure. This choice is indeed correct when only minimal 
sensors are taken into consideration. M. Uchiyama and K. Hakomori proposed the use of 
the condition number of the compliance matrix of a force sensor in [19]. Unfortunately, their 
work was published in Japanese, and it has not been as widely known as it deserved until 
cited by [20] 2 . The authors of [20] criticize the condition number criterion of [19] under four 
regards. Some of those points are discussed here, since this is believed to provide insight in 
the general problem of optimal design of force sensors. 

The first remark is that the condition number criterion is insensitive to absolute values 
of the compliance matrix entries. In other words, N( C) = N(2C). While it is an obvious 
advantage of condition numbers to be scale-independent, the authors of [20] underscore the 

2 An independently developed formulation of the method has been presented at the same time in [11] 
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fact that sensors with larger absolute values of structural strains are not rewarded by this 
criterion. A design principle consisting in making the “strain-gauge sensitivity” as large as 
possible is proposed, which consists in maximizing the norm of the rows of C. In other words, 
a sensor is optimal under this principle if, for the maximum nominal load, every gauge is 
strained at the maximum allowable level. 

The second, third and fourth remarks to the condition number criterion all apply to the 
fact that it does not generalize to redundant sensors. For instance, the two compliance ma¬ 


trices Ci = 


1 0 
0 1 


and C 2 


1 0 \ 
1 0 

0 1 


are brought to our consideration. The condition 


numbers are Nq 1 = 1 and Nc 2 — \/2, respectively. The authors note that, although C 2 may 
have a redundant strain-gauge, C 2 seems not to have less sensing performance than Ci . 
This fact is true indeed, and it will be shown in the following sections that this problem 
(as well as the others pointed out by [20]) can be solved by a correct generalization of the 
condition number criterion to redundant sensors. 

Two more design principles are proposed by the authors of [20], the first being to maximize 
the minimum singular value of the compliance matrix C, and the second to minimize the 
maximum singular value of a displacement matrix G, evaluated at the point of application 
of the load. 

In the method presented in this paper we avoid the definition of multiple optimality 
criteria, in order to avoid the need of a difficult choice of weights, which always imply a degree 
of arbitrariness. A lower bound on sensitivity, as well as an upper bound on displacement 
under load, can be very easily incorporated as unilateral constraints to the design (see 
section 2 above and 7 below). On the other hand, an attempt is made at generalizing the 
condition number criterion for redundant sensors, as the single criterion able to synthesize 
accuracy requirements for the sensor (it can be seen for instance that maximizing strain 
gauge and force sensitivities as defined by [20] leads to minimizing the condition number). 


5 Redundant Sensors 

Consider again equation 1, v = Cp, in the case that the number of independent measure¬ 
ments, m, is larger than the number of unknown load components, n. Since the entries of 
both C and v derive from experimental measurements, equation 1 is in any practical case 
inconsistent, i.e. no exact solution is possible. Only approximate solutions, denoted with p, 
can be obtained through a manipulation of equation 1 as: 


(5) 

( 6 ) 
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Mv = MCp, 

p = (MC)" 1 Mv, 




Figure 1: Error propagation in a force/torque sensor 

where M € 3? nXm , and it is assumed that det(MC) ^ 0. 

Such a manipulation consists of projecting the measured point v of the m-dimensional 
space onto the range space of C, denoted with 7£(C), which is a n-dimensional subspace. The 
direction along which the projection is made is specified by the choice of M. In particular, 
if M = C T , the projection results normal to 'R.(C), and the corresponding solution p = 
C^v (where ct = (C r C) -1 C T is the pseudo-inverse of C) is the least squares solution 
to equation 1. Note that we used here the assumption that C is full rank, which is not 
restrictive at all insofar as reasonable sensor designs are considered. 

Naturally, any minimal sensor can be considered as a particular case of a redundant one, 
by supposing that (m — n) measurements are neglected. Correspondingly, and modulo a 
reordering of rows, a manipulation matrix M = (Inxn|0nx(m-n)) is adopted in equation 5. 
Also for minimal sensors, therefore, the solution p is to be considered only an approximation 
of the “true” value of applied load. 

The total error £p can hence be thought as a sum of two terms, deriving from the 
projective manipulation 5, and from the actual solution of 6, respectively. The rest of 
this paper will discuss the application of the generalized condition number criterion to the 
solution of equation 6: minimizing the error in its solution is assumed as the goal of design 
optimization. Before beginning such discussion, however, it seems appropriate to show that 
the error term deriving from the projection 5 is not directly affected by the choice of the 
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C) 


Figure 2: The projection error is not affected by the design of the compliance matrix. 


compliance matrix C elements, but rather it depends only on the source errors, Sv and 
SC. In other words, the projection error is technology-dominated, while the solution error 
is design-dominated. 

This fact can be graphically illustrated with reference to a very simple example with 
m — 2, n = 1. An hypothetical error-free sensor, modeled by the linear equation v = Cp, 
would have the measurement vector v lying exactly on 72-(C), represented in figure 2 with 
a dashed line. The “true” solution would be represented in this case by the length of the 
segment OZ. However, any solution of the real sensor sketched in equation 2 corresponds to a 
projection of the vector v Sv onto 7?.(C + <5>C) (solid line in figure 2). Different manipulation 
matrices produce different approximate solutions. For instance, the length OA corresponds 
to the choice M = (1 0), i.e. to neglecting the second measurement; OC corresponds to 
M = (0 1), neglecting the first one; OB represents the least-squares approximation obtained 
with M = C T . The “true” solution, of course unknowable, can now be represented by the 
length OY. The total error 6p is the sum of a projection error (e.g., OB - OY in the case 
of least squares), and of an algorithmic error, which propagates in solving equation 5 (i.e. in 
finding the length OB). 

The sensor design, i.e. the choice of the elements of C, place the subspace 7L{C) in 
3? 2 . By varying this choice, and keeping Sv and SC constant, the projection error does not 
change, since this variation amounts to a rigid rotation of the graph of figure 2 about the 
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point 0. Therefore, it seems reasonable to apply design criteria that tend to optimize the 
accuracy in evaluating p , the best affordable approximation of the “true” load with a given 
technological implementation. 

Applying equation 3 to equation 5, we obtain: 


€p — ^mc)-A"p(hdC), 


(7) 


where e r 


ll<(Mv) 

||Mv|| 


^mc 


||5(MC)|[ 

||MC|| 


In order to correctly compare sensors that employ different solution methods for equa¬ 
tion 5, or even sensors built with different numbers of basic transducers but comparable 
technology, we need to express equation 7 explicitly in terms of relative source errors e v , 
e c . By “comparable technology” it is meant here that all basic transducers have the same 
statistical properties in terms of measurement and calibration accuracy, such that e v and e c 
are not varied by adding or removing transducers 3 . 

The term e mc can be rearranged as 


€mc 


||(M + £M)(C + £C)-MC|| 
l|MC|| 

|[M <j>C + <j)M C|| < ||M|| |1C 
||MC|| _ ||MC|| 


0, 


where e m = is the relative error the matrix M is known with. As for the term e mv , it 
m ||M|| 

holds: 


e 


mv 


||(M + 6M)(v + 8v) — Mv| 
l|Mv|| 

||M 6v + SM v11 ^ ||M|| ||v 
||Mv|| “ ||Mv|| 


~ ( e m + tv)- 


( 8 ) 


In order to give a proper upper limit to the last term of equation 8, which propagates 
source errors into the computed measurements Mv, it is expedient to consider only vector 
and matrix 2-norms (i.e. the Euclidean and maximum singular value norm, respectively). 
It should be pointed out that, although this choice was not specified before for the sake of 
generality, the use of 2-norms is always advisable in problems like the present one. In fact 

3 There may actually be an improvement of by adding more transducers with the same individual 
accuracy, in the sense that more measurements of the same quantities would be available. The same effect 
could be obtained on the other hand by repeating the measurement of some or all transducers and averaging 
before solving the equations. We assume for brevity’s sake that both those effects are accounted for in 
evaluating e„. 
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these norms have the property of being invariant with respect to orthogonal transformations, 
due e.g. to rigid rotations of the reference frame. On the other hand, being the dimensions 
of vectors and matrices involved in force sensing small, the computational load is not much 
larger for 2-norms than for other possible norms. In the following only 2-norms will be 
considered, and explicit notation will be omitted. 

The ratio 


+ 


< 


IMI 


l|Mv|| 


< 


IMI 


infi 


IlMvIl 


(9) 


mi ght, appear unbound at varying v. In particular, for v lying in the null space of M, A/"(M), 
infinite amplification of source errors would occur. This corresponds to a choice of M that 
neglects every meaningful information from the sensor. More precisely, an upper bound to 
the ratio 9 can be obtained as follows. Consider the orthogonal decomposition of the vector 
v = Vp+Vfc, with v h G M{ M), and v p <E K{M r ). Using equation 1, we can write = C p*. 
This would not be possible if v/, ^ 7£(C): but in such case V/, would not correspond to any 
physically meaningful load configuration, so that it must be € TZ(C). On the other hand, 
if we apply equation 5 we have that 


Mvfc = 0 = MCph. 

This equation, along with the previous assumption det(MC) ^ 0, leads to the conclusion 
p h = v h = 0. 

Since for what has been shown the existence of components of v belonging to Jv (M) can 
be excluded, equation 9 can be modified to 

c™ < 1|M|| = ||M|| 

Cm + € v infvj, Pmin 

where /i min indicates the smallest non-zero singular value of M (i.e., the square root of the 
smallest non-zero eigenvalue of M r M). It can be easily verified that 

ML = IIMM^II = N( M) = JV M , 

ft min 

where the obvious generalization of the definition of 2-norm-based condition number to 
rectangular matrices is introduced [21]. 

Finally, equation 7 can be restated as 

c p < + £m) + N M (e v + c m )| K p { MC) = F(C,M,e v ,e c ,e m ). (10) 
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This formula is the most general result of this paper. The design objective function F(-) 
is defined as a function of the compliance matrix C, of the adopted manipulation M, and 
of the relative source errors e c , e v , and e m . Note that C and M are in turn functions of the 
design variables x. Most often Cm can be expressed in terms of C and e c , so that it will be 
dropped from the arguments of F(-). Given the values of relative source errors e v , e c that 
are expected with the adopted degree of technological sophistication, the accuracy of force 
sensors differing in design, solution method and number of transducers (provided that all the 
transducers have similar statistical properties) can be compared by means of equation 10. 

Two particular cases are very important in practice, and are briefly discussed in the 
following subsections. 

5.1 Sensors Solved With Least Squares Methods 

A general least-squares approximation of the inconsistent set of equations of a redundant 
force sensor can be obtained if the choice M = is made in equation 5, where is 

a diagonal matrix of weights corresponding to different accuracies of basic measurements. 
Equation 10 can be applied to estimate error propagation in such a sensor. However, since 
a correct weighting is rarely achievable, the regular least-squares solution with M = C T is 
the only one currently adopted in practice (see e.g. [15]). 

Applying equation 10 and considering that, for 2-norms, it holds ||C r || = |]Cj| and 
||C T C|| = ||C|| 2 , we have: 


€ P < {[2 + N c ] e c + N c e v } K p (C T C), (11) 

where 


K p { C T C) = 


Nl 


1 


2 N 2 c e c 


Nc- 


( 12 ) 


Equation 11 and equation 12 show how the dependence of sensor accuracy on the condi¬ 
tion number is increased for sensors using least-squares approximation 4 . On the other hand, 
redundant transducers can reduce the condition number Nc- In summary, the effectiveness 
of redundant transducers should be checked in the specific case, as they may result either 


4 It can be noted that equation 11 and equation 12 suggest a more pessimistic estimate of error propagation 
in least square problems than is currently held in numerical stability analysis [21]. Our bound on error 
propagation is in effect quite conservative; propagation problems with redundant sensors can be relaxed using 
more sophisticated algorithms for the inversion of the normal equation 5 with M = C T , e.g. Householder 
or SVD methods [23]. More precise bounds would however involve an a posteriori evaluation of errors, 
which is not clearly feasible in our problem. Relations 11 and 12 provide a concise objective function for 
optimization procedures, whose degree of conservatism is probably of the same order for different designs 
put in comparison. 
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beneficial, or useless, or even harmful for sensor precision (this point will be reconsidered 
later). 


5.2 Sensors With Structured Matrices 

Some force sensor designs produce compliance matrices whose elements obey to such a pat¬ 
tern that the solution of the associated linear system can be performed very easily. Well 
contrived sensor designs in fact can realize diagonal or near diagonal MC matrices, thus de¬ 
coupling the effect of each load component on the manipulated measurement vector Mv. A 
good example is described in [5]. Although this concept is not limited to redundant sensors, 
it is easier to design decoupled sensors when m > n. 

Good decoupling was a strict requirement for force sensors designed before the advent 
of present cheap and powerful electronic computers. In fact, all the data manipulations 
needed for such sensors could be carried out by simple analog summers and multipliers, 
most often realized by arranging the Wheatstone bridges and the amplifier gains used for 
measurement in a peculiar combination. In some cases, for instance when requirements 
on sensor bandwidth dominate over accuracy needs, such designs can be still today used 
effectively [20]. 

The M manipulation matrix for decoupled sensors consists of a simple pattern of integer 
numbers, so that we can assume ~ 0. Equation 10 can therefore be simplified to the 
form: 


+ (13) 

Some interesting conclusions can be drawn from equation 13. It can be noted for instance 
that a rescaling of the sensor equations, tending to equilibrate them by multiplying the 
characteristic matrix by a diagonal weighting matrix M (which obtains a lower condition 
number for MC), does not produce any enhancement of sensor accuracy. This fact has 
obvious physical meaning, and is in effect equivalent to Bauer’s theorem, well known in 
numerical analysis 5 [22]. 

Furthermore, if sensors are considered as particular redundant sensors, we expect that 
the objective function calculated with equation 3 should agree with the one resulting from 
equation 10, in its specialized form 13. This can be in fact easily verified by substituting 
M = (lnxn|0„x(m-n)) in equation 13. 

Finally, equation 13 allows us to address the problem raised in section 4, i.e. that the two 


matrices C x 


1 0 
0 1 


and C 2 


1 ° \ 
1 0 

0 1 / 


have different condition numbers, although 


5 Note that this scaling is purely algebraic, as opposed to a weighting based on statistical properties of 
measurements, which could in fact increase the sensor accuracy, affecting the source error e v 
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the same information seems, so to say, to be stored in them. The point, which has general 
validity, is that it makes no sense to consider the objective function value associated with the 
sensor C 2 , without specifying how the equations are going to be solved, i.e. which matrix M 
is to be used. So, if we decide to use least squares methods, application of equation 11 with 
e c = 0 gives the objective function value .F(C 2 , C^, e„, e c ) = 2\/2e„, which is much worse 

than F( C X ,I, e„,e c ) = e v . Yet, a manipulation M = ^ J jj J ^ can be applied, such that 

-F(C 2 , M, e v , e c ) = e„. 

Among the factors that tend to make the decoupling approach obsolete in force sensor 
design, one is the complexity of mechanical structures necessary to implement it. Design 
complexity conflicts with requirements on size, cost, and reliability of the sensor. Further¬ 
more, and more subtly, the complexity of structures (where for instance stress concentrations 
are more or less deliberately introduced), makes the linear model on which the sensor is based 
even more arbitrary. Besides that, it should be pointed out that force sensors have structured 
compliance matrices only as far as they are on the designer drawing table. No calibration 
technique is allowed by decoupling design to correct implementation dependent errors on C, 
so that larger errors e c must be expected (typically an order of magnitude larger than in 
calibrated sensors, in the author’s experience). 


6 Accuracy And Adding More Transducers 

Equation 10 identifies an objective function capable of leading an automatic design proce¬ 
dure in the choice of optimal solutions. Nevertheless, a better understanding of transducer 
redundancy can be gained by analyzing how Nc is varied by adding new transducers. 

If C € 5? mXn is the compliance matrix of a given (stage in the design of a) sensor, 
an additional set of d basic transducers (henceforth called gauges) will produce the new 

characteristic matrix A = ^ ^ , where D € 3i <ixn , A £ sft( m+<i ) xn . The condition number 

of the new linear system, according to the choice of 2-norms, is the square root of the ratio 
between the largest and the smallest eigenvalues of A r A. Considering A . T A as the sum of 
two symmetric matrices, 

A r A = ( C T I D 3, ) f-g-) = C T C + D r D, 

and applying the minimax theorem [23], we have: 

max t - = i, n { 7 ? + < N 2 < 7 1 + S i ( 14 ) 

min i=1>n { 7 ? + 6„_ i+ i} ~ A ~ ll + ’ 
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where ( 71 , 72 , • • • , 7 n) and (8 1 ,8 2 , ...,8 n ) are the singular values of C and D in non increasing 
order, respectively. This relationship gives bounds for condition number variations as a 
consequence of the addition of new gauges, but it is insufficient to determine whether the 
variations are favorable or not. 

More detailed information can be obtained if extra gauges are added one at a time. The 
matrix D in this case is reduced to a row-matrix d r £ 3? lxn , and D r D = d r d £ 3ft nXn has 
singular values £1 = (d^d) 1 ^ 2 , S 2 = £3 = • • • = 8 n — 0. Any set of n orthogonal vectors which 
includes d forms an eigenvector set for D r D. Equation 14 is modified in this case as 

max{7 2 ,7 2 + ff} ^2 < 7 i+_£i 
m i n {7n-n7n + *l} “ A ~ ll + 8 l 

If d is chosen parallel to the eigenvector of C T C corresponding to the smallest eigenvalue 
Ini C T C and D r D share the same eigenvector set. Therefore, also A T A = C T C + D T D 
shares the same eigenvectors, and the eigenvalues of A T A are equal to those of C r C except 
for one, a] = 7 2 + d T d, with 1 < s < n. 

Corresponding to this choice of the direction of d, we have the new condition number: 


maxi =1>n {7^ 5 7 n + dTd } 
min i= i >n {7n—7n + d T d} 


According to this result, a new transducer should be placed on the sensor in a position 
such that its output v m +i is related to the load p as v m +i = d r p, with d parallel to the 
singular vector corresponding to the minimum singular value of the old characteristic matrix, 
and such that d r d = 7? — 7„- By satisfying perfectly these conditions with n — 1 additional 
transducers, it is theoretically possible to modify any given design so as to obtain unitary 
condition number. In practice, this goal may be impossible to achieve, due to the fact that an 
arbitrary output law can not always be obtained from a transducer. However, these design 
suggestions, along with equation 10, can be useful to decide whether, and how, to correct 
sensor designs by adding transducers. 


T Applications 

In order to illustrate some aspects of the methods above outlined, two applications of op¬ 
timal force sensor design are presented in this section. The first one is extremely simple, 
and intended only to illustrative purposes. The second example reflects a real application 
problem, regarding the design of miniaturized force sensors designed to fit robotic fingertips. 
The sensor is currently being produced and successfully employed in several applications. 
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Figure 3: A simple cantilever beam sensor with rectangular cross section. 

7.i A Simple Two Axis Force Sensor 

An extremely simple two axis force sensor can be realized as depicted in figure 3. By means 
of strain measurements on the cantilever beam with rectangular cross section b x h and 
length /, the unknown load components p\ and P 2 are sensed. The five components of the 
vector x of design variables are defined as Xi = b, = h, and X3,X4,Xs as the distances 
of the gauges from the tip of the beam. Three unilateral constraints express limitations on 
acceptable stress and displacement under the maximum nominal load p max , as well as on 
the strain under the minimum nominal load p m ; n : 

s,(x) = + £*==51 _ ^ < o, 

Sj( *) = (16) 

/ X Pi ,mtn . Pi ,mtn 

03(x) =- - -1- - — -h Un < 0, 

where A = is the cross section area, J = its momentum of inertia, E is the 

Young’s modulus of elasticity of the material, and cr adm its maximum admissible strain; 
fmax is the maximum acceptable displacement under load, and £ mtn the minimum sensitivity 
required to the sensor. The feasibility region is defined as 
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0 < Xi < l, i = 3,5, 

while obvious bounds are set on ij = i and x 2 = h by constructive considerations. Three 
cases will be considered: a minimal sensor solved with full inversion or with special solution, 
and a redundant sensor with three basic transducers. 

7.1.1 Minimal realization with full inversion 

For a minimal sensor using only two gauges, Si and s 2 , the strains are: 

Pi t max Pl . P2,max x 2 P2 

Vl = - 1 -—- *3 -, 

EA P\ t rnax EJ P2,max 

Pl ,moi Pl P2 ,7Tl <IX x 2 P2 

V 2 = ——-—— x 4 -• 

&A. Pl,max p 2 >max 

The sensor model is v = Cp, where p = '-— —^ —) is the normalized load vector, 

\Pl,max P2,max / 

and the compliance matrix (modulo a constant) is 

( Pl ,max P2 t max &2 ^3 \ 

EA EJ ) 

Pl ,max P2 ,max'^2 J * 

EA ~ EJ / 

In this simple case, the optimum design can be easily found. Imposing Nq = 1 amounts 
to an orthogonality condition for the columns of C, c x and c 2 : 


whence the design 


cfc 2 = 0, 

ll C l|| = ll C 2 


*^3 — — 


Pl t max 
12j2 P 2 ,max 


Given such dimensions the resulting error is e p ps e v + e c . Constrained minimization 
techniques should be used if these dimensions do not comply with constraints 16. 
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7.1.2 Minimal realization exploiting the structure of C 

If the peculiar structure that C has for X 3 — = x is exploited by using a manipulation 

matrix 


M = 


1 1 

1 -1 


the original equation is modified to 

' n + v 2 A 


Mv 


0 


Vi - v 2 


( p Pl, max 

q A 2 P2 - maxX2X 


which can be solved with only n = 2 flops, instead of the n 2 = 4 otherwise necessary. 
This advantage is obviously of some, if little, importance only for larger n than used in this 
example. 

The following relationships can be easily verified: 


l|M|| = ||M 


-H 


■\/2; Nm = i; 

,m<*x x 2® 1 


I^ill P f Pi,max P2,t..— _ . . 

| C || = vWj—, — 


Itk /rII o f Pi,max Pl,max x 2 x \ 

| M C|| = 2max(— , —gj-}, 

Jpi 

{max A-P 2 ,m ax 2 ® 


N{ MC) = max 


Ap2 t max*£2*£ Jpl^n 


According to equation 13, if x = is chosen, we have 

' IMII IICI 


e p < 


IMCI 


-e c -f Nm^v f ifp(MC) « e v + e c , 


which shows that this manipulation does not worsen error propagation. However, as al¬ 
ready mentioned, in practical cases the cancellation of the off-diagonal entries of MC is not 
expected to be perfect, so that larger e c must be taken into account. 

7.1.3 Redundant realization 

If the three gauges si, s 2 , s 3 are placed on the beam, the redundant sensor compliance matrix 
results 



/ Pi,max 

P2,max*2*3 


EA 

EJ 


Pi ,max 

P2,max x 2 x 4 


EA 

EJ 


Pi ,max 

P2 ,max«^2 *^5 


\ EA 

EJ 
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By imposing Nc = 1, we obtain the optimal design: 


T 

2 i 2 i 2 

*3 d - * 4 + a : 5 

and, applying equation 11, we have 

Cp ^ “i“ ' 

This result shows that the introduction of the third gauge makes error propagation worse. 
On the other hand, the optimality conditions above allow a wider design flexibility, which 
may be necessary, depending upon the constraints applied to the variable ranges. 

7.2 A Six Axis Miniaturized Force Sensor 

As mentioned in the introduction, a particularly promising application of force sensors is 
intrinsic contact sensing [8]. In its most common embodiment, this method adopts six-axis 
force sensors integrated as close as possible to the parts of the arm whose surfaces contact 
the environment. This application has critical size requirements: for instance, force sensors 
have been designed to fit the interior of fingertips of articulated dextrous hands as small as a 
cylinder 18mm diameter and 25mm long. Furthermore, contact sensing algorithms use force 
sensors outputs as input data, thus making accuracy requirements also very important. 

The approach to the design of sensors for intrinsic contact sensing that will be illustrated 
here basically follows the guidelines illustrated below, which we consider in general advisable 
for high-performance sensor design: 

• Use very simple mechanical structures, whose behavior is as close to linear as possible; 

• Utilize the least necessary number of transducers, to avoid unnecessary propagation of 
errors; 

• Give the design a convenient parameterization; 

• Use numerical techniques to find the best combination of design variables which com¬ 
plies with design constraints. 

Following this approach, a miniaturized force sensor has been designed [24] and built 
in several copies, used to sensorize devices for dextrous manipulation such as the Salisbury 
hand [25], the “tactile explorator” finger of the Centro “E.Piaggio” [11], and the whole-hand 
manipulation system UB-Hand II [3]. The mechanical structure of this sensor is sketched 


= 0 

3 Pjmax^l 
144 pl max 
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Figure 4: A miniaturized force/torque sensor with cylindrical beam structure. 

in figure 4: it consists of an hollow, thin-walled cylinder. Strain-gauges are applied on the 
external surface of the cylinder, in a number of six 6 . The cylinder dimensions, the position 
of the gauges on its structure and their orientation have been considered as the variables to 
be optimized. 

Figure 4 shows the {0, Ci,C2>C3} reference frame in which the components of the load 
applied to the sensor extremities are described. The axis (i is placed along the cylinder axis. 
The position of the i th gauge is uniquely determined by the cylindrical coordinates of its 
center point, and 0,, and by the angle <f> t formed by the gauge axis with the cylinder axis. 
The design variables are three for each gauge, plus the cylinder radius and wall thickness, 
i.e. 20. 

The simplicity of the mechanical structure allows the evaluation of the entries of the 
compliance matrix C by means of simple relations of elastic beam theory [26]. Note that the 
computation of C has to be done at each iteration step of the numerical optimizing routine, 
so that avoiding demanding techniques as the finite element method is extremely expedient. 
Accordingly, the strain corresponding to the i th gauge is 

v i ~ y GujPji 

j= 1,6 

6 This number is actually increased to 7 for temperature compensation; this point is addressed in the 
appendix. 
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where Pj is the value of the j th load component normalized with respect to its nominal 
maximum value, Pj tTnax . Elements are as follows: 


Ci,i 

— Pl,max 

W n (cos 2 4>i~ v sin 2 

4>&, 



C it2 

— P2,max 

[Wf(i cos ^(cos 2 <f>, 

i — i/sin 2 <f>i) + 

W, sin 6i 

sin2^j ; 

Ci,3 

— jP3,maa? 

[WfQ sin 6i( cos 2 4>i 

— i/sin 2 <j>i) + 

W 3 cos 0 % 

sin2<^ij ; 

c iA 

— P^max 

W t sin 2 <j>i\ 




Ci,s 

— P 5 tin ax 

Wf cos ^j(cos 2 <f>i — 

i/sin 2 (f>i ); 



C',6 

~ P6,max 

Wf sin 0;(cos 2 <f>i — 

i/sin 2 fa); 




where v is the Poisson’s ratio for the structure material, and the moduli W are defined as 


1 

2irRsE ’ 

2Wn . 

R * 

2(1 + u)W n ; 

(1 + u)W n 
R 

where E is Young’s modulus, R is the cylinder radius and s its wall thickness. The relations 
above are valid in the assumption s <C R: in this case, l/s is a common factor of Cij, so 
that Nc is not affected by s. In other words, the relative accuracy of the sensor does not 
depend on the wall thickness; s can be chosen independently of other design variables, and 
made such that minimum sensitivity and maximum strain level requirements are met. 
Other design constraints for a sensor to fit a robotic hand fingertip are of the form 


W n = 

Wf = 
w, = 
w t = 


o > \Ci~Ci\~ ^ TflCtX ) 

Rmax > R>2S, 

where l max = 15mm, Rmax — 6mm are typical values. 

According to the previous discussion, an objective function F(R., <f>i) = Nc has been 

assumed to optimize this minimal sensor. Note that for minimal sensors, we do not need to 
use explicit values of the source errors e„ and e c . The condition number of the compliance 
matrix is a complex function of the 19 design variables whose minima could not be found 
analytically. A numerical algorithm (a customized version of Powell’s and Broyden’s methods 
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[22]) has been employed. The values of design variables corresponding to the lowest of the 
relative minima found in repeated trials (resulting in Nc, m in = 3.07), are given in table 1. 
A singular value decomposition of the corresponding compliance matrix is reported in table 
2, showing that a 3.07 ratio occurs between the maximum and minimum relative sensitivity, 
the latter corresponding approximately to the pi component of the load, acting along the 
cylinder axis. 


Variable 
Gage no. 

</>(rad) 

0(rad) 

C(mm) 

iZ(mm) 

s (mm) 

1 

0.3 

-0.7 

1.2 



2 

-0.3 

1.6 

1.7 



3 

0.2 

3.2 

-1.9 

5.4 

0.2 

4 

-0.6 

-1.8 

-1.4 



5 

0.4 

0.3 

-5.9 



6 

0.3 

0.9 

6.0 




Table 1: Optimal design variables for the six-axis miniaturized force sensor, corresponding 
to the nominal load: Pi = p 2 = Pz = 10JV; p± = 607Vmra; ps = Pe — 85 Nmm. 



Singular Structure 

Singular 

7i 

72 

73 

74 

75 

76 

values 

138.4 

129.6 

89.4 

64.3 

55.9 

45.0 

Singular 

Ul 

u 2 

U 3 

u 4 

u 5 

u 6 

vectors 

0.07 

0.08 

0.00 

-0.05 

-0.27 

0.95 


-0.04 

0.32 

0.83 

-0.03 

-.17 

-0.04 


0.05 

-0.36 

0.10 

0.81 

-0.44 

-0.05 


0.08 

-0.18 

0.31 

0.34 

0.82 

0.26 


0.02 

0.83 

-0.28 

0.46 

0.10 

-0.02 


0.91 

0.16 

0.34 

-0.10 

-0.11 

-0.12 


Table 2: Singular values and corresponding singular vectors of the sensor compliance matrix 
for the optimal set of design variables listed in table 1. 
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Figure 5: Sub-optimal arrangement of gauges on the miniaturized cylindrical sensor. 
7.2.1 Sub-optimal design 

An interesting sub-optimal configuration of the same force sensor is depicted in figure 5. The 
centers of the gauges are placed on two circles, on the plane £i = 0 for gauges 1 through 
4, and on the plane = —d for gauges 5 and 6. Angles are $i = z|, and <f>i — ( —1 ) l 4>, 
1 < z < 6. Design variables are thus reduced to three, R, d, and^>. The compliance matrix 
given by such design has the following structure: 

/ a 0 -b e 0 / \ 

a b 0 —e / 0 

_ _ a 0 b e 0 —/ 

a —b 0 — e — / 0 

a c —b e 0 / 

\ a b -c -e f 0 , 

An almost-decoupled sensor can then be obtained even with this extremely simple design, 
by using a manipulation matrix 
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/ 

M = 

\ 


1 1 

-1 0 

0 1 

1 -1 

0 1 

1 0 


1 1 

0 0 

0 0 

1 -1 

0 -1 

-1 0 


0 0 \ 
1 0 
0 -1 
0 0 
0 0 
0 0 / 


An extremely fast solution of the sensor equations is possible with this arrangement: 


f Vl + v 2 

+ V 3 

+ v 4 \ 


f 4a 

0 

0 

0 

0 

0 

\ 

v 5 - Vi 




0 

c 

0 

0 

0 

0 


V 6 — V 2 




0 

0 

c 

0 

0 

0 


Vi - v 2 

+ V 3 

- v 4 


0 

0 

0 

4e 

0 

0 


V 2 - v 4 




0 

2 b 

0 

0 

2/ 

0 


t— * 

1 


) 

1 1 

l 0 

0 

-2b 

0 

0 

2/ 

) 


The optimal combination of design variables with sub-optimality constraints is reported 
in table 3. The corresponding objective function minimum is at Nc = 6.6. As expected, the 
sub-optimal sensor is less accurate than the optimal one. The advantage of this design is 
that it is possible to choose either a precise solution algorithm (full inversion of C) or a fast 
one. Perhaps more importantly from a practical point of view, the symmetries of the design 
render its fabrication much easier than the fully optimized one. It is finally noted that the 
objective function is fairly smooth in the neighborhood of its minimum, which fact renders 
the design quite robust to fabrication inaccuracies. 


Variable 
Gage no. 

^(rad) 

0(rad) 

C(mm) 

J2(mm) 

s(mm) 

1 

0.3 

0 

-7.9 



2 

-0.3 

7v/2 

-7.9 



3 

0.3 

7T 

-7.9 

5.6 

0.2 

4 

-0.3 

3tt/2 

-7.9 



5 

0.3 

0 

4.1 



6 

-0.3 

7r/2 

4.1 




Table 3: Sub-optimal design variables for the six-axis miniaturized force sensor, correspond¬ 
ing to the nominal load: Pi = p 2 — p 3 = 10A^; p 4 = 60 Nmm\ p 5 = p& = 85Nmm. 
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Singular Structure 

Singular 

7i 

72 

73 

74 

75 

76 

values 

121.4 

121.2 

65.0 

60.2 

32.6 

32.2 

Singular 

Ul 

u 2 

U3 

u 4 

u 5 

u 6 

vectors 

0.00 

0.07 

-0.54 

0.00 

0.83 

0.00 


0.42 

0.42 

0.49 

-0.51 

0.28 

0.25 


-0.42 

0.42 

0.49 

0.50 

0.28 

-0.25 


-0.18 

0.00 

0.00 

0.31 

0.00 

0.93 


0.55 

0.56 

-0.34 

0.44 

-0.26 

-0.04 


-0.55 

0.56 

-0.33 

-0.43 

-0.26 

0.04 


Table 4: Singular values and corresponding singular vectors of the sensor compliance matrix 
for the sub-optimal set of design variables listed in table 3. 


7.2.2 Some Experimental Data 

As it has been pointed out several times already, the ultimate performance of a force sensor 
depends on design quality as well as on technological factors. It is not therefore easy to assess 
design quality by experimental means, as an extremely large statistical basis, and a tight 
control of experimental conditions (including details such as quality of gauges, of bonding 
agents, of electronic components and so on), would be required. Although such exhaustive 
verification can not be provided here, some numerical values obtained from miniaturized 
multi-axis sensors will be given for generic reference. 

An aluminum (2024 T4) sensor designed according to the sub-optimal scheme, with wall 
thickness 0.8mm and using foil-gauge components, has been connected to instrumentation 
amplifiers and to an A/D port of a personal computer. The maximum nominal load was 
30N for force components and 150Nmm for moment components. An estimate of the relative 
source errors is to ~ e c ~ 1%. As a figure of the global accuracy of a sensor, the cross-talk 
matrix X is sometimes used, whose entry is defined as the ratio between the measured 
value of the i th component of the nominal load, p;, and the actual value pj (for an ideal 
sensor, X = I). The cross-talk matrix obtained for the sensor (employing a full inversion 
solution scheme) is as follows: 
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/ 1.02 

-0.01 

0.001 

-0.01 

0.01 

0.01 

\ 

0.01 

1.01 

-0.03 

-0.01 

0.00 

0.00 


0.01 

0.01 

0.99 

0.05 

0.00 

0.01 


0.01 

-0.02 

0.00 

0.98 

-0.01 

-0.02 


0.00 

-0.01 

-0.02 

-0.03 

1.02 

0.01 


\ -o.oi 

-0.01 

-0.01 

0.04 

0.04 

0.99 

1 


Assuming the maximum norm of the columns of X as a measure of the relative error, 
e p = 7.5% can be obtained by this matrix, which is lower than the expected bound e p < 
N c e v e c ~ 13%. This result is believed to depend partly on the conservative nature of the 
condition number criterion, partly on the fact that the “worst possible combination” of load 
components is not reflected in the cross-talk matrix. 


8 Conclusions 

It has been shown that the design of a multi-axis force sensor can be optimized with respect 
to its accuracy by minimizing an objective function of the design variables. The proposed 
function of merit (equation 10) is general enough to allow the choice of the optimal design in 
a broad range of possible solutions. In particular, the choice of the optimal number of basic 
transducers to be used in a multi-component sensor is addressed. For sensors employing as 
many basic transducers as the unknown load components (minimal sensors), the objective 
function reduces to the condition number of the sensor compliance matrix (equation 3). 

An often asked question about multicomponent sensing is: “How many basic transducers 
should be utilized? Are redundant transducers useful, or useless, or are they harmful to the 
overal sensor performance?” A clear-cut answer to this question is not possible, since in the 
general case the direct evaluation of equation 10 is necessary to assess the results. However, 
a general guideline that results from the discussion of equation 10 is that, whenever the 
optimal configuration of a minimal sensor complies with the constraints, adding redundant 
sensors is unlikely to improve sensor accuracy. 

Finally, we would like to note that, although multi-axis force sensors have been explic¬ 
itly considered in this paper, the method proposed is more general, and can be directly 
extended to any sensor of multiple physical quantities based on a number of linearly related 
measurements. 
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Appendix: Compensation of Thermal Effects 

The various effects of temperature variations are among the most important error sources 
for basic strain transducers. For resistive strain-gauges, such effects can be grouped, in 
decreasing order of importance, as follows: 

• variation of the intrinsic (non stressed) gauge resistance; 

• thermal expansion of the sensor structure; 

• gauge-factor variation; 

• variation of the elasticity modulus of the structure material. 

While the first two phenomena cause errors in measurements v;, the latter two modify 
the compliance matrix elements, entering the measurement output in a nonlinear manner. 
Strain gauges exhibiting a curve of intrinsic resistance vs. temperature that approximately 
matches (at least in some temperature range) the thermal expansion vs. temperature char¬ 
acteristic curve of the material, are currently available. Such gauges perform a first, rough 
compensation of thermal effects. 

Possible approaches to finer compensation are listed as follows, in increasing order of 
accuracy: 

• To regard the effects of temperature variations AT as a component of the source errors. 
The counteraction is to minimize error propagation; 

• To regard nonlinear effects of temperature as source errors, an assume that the remain¬ 
ing effects are functions of the temperature only, independent of the individual gauge 
characteristics, position, orientation etc. Since the effects of temperature are common 
to every gauge on the sensor, a dummy gauge applied on a stiff part of the structure, 
and subject to the same temperature variations, can be taken as a reference. 

• To consider nonlinear effects as source errors, and model the remaining effects of tem¬ 
perature T with suitable individual functions r,(T). These functions are typically 
approximated by polynomials of degree three or higher. Starting from data normally 
provided by the gauge supplier, the width Wb(T ) of a band within which the r; of every 
gauge in the batch employed lies, can be obtained (see figure 6.a). If two curves r/(T) 
and r u (T ) are considered such that, for any T and any i , r/(T) < r,(T) < r u (T), we 
have: 
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Figure 6: Typical shape of foil gauge resistance vs. temperature curves: a) absolute; b) 
relative to an offset measurement. 


r u (T) - ri(T ) < w b r u (T). 

If moreover an offset measurement is taken at T = T 0 , the curves Ar^(T — T 0 ) = 
Ar^AT), can be represented as in figure 6.b. Let Ar m (AT) = Ar >+ Ar » be the median 
curve of the band, so that any curve Ar, can be approximated by Af; ~ ct t iAr m (ST), 
where c tj ; is a constant characteristic of the i th gauge. An upper bound for the approx¬ 
imation error is w b . The effects of the temperature on a single measurement can be 
expressed as 

v { = CiP + c t)l Ar m (AT), 

where C; is the i th row of C. Equation 10 can therefore be rewritten as 
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The elements c (it of the column vector c t can be calibrated in controlled temperature 
conditions. The compensated sensor requires the solution of a linear system of m 
equations in 7 unknowns. The seventh unknown, Ar m , is actually of little practical 
interest, and its calculation can be superseded, as long as the whole calibration matrix 
has been used to evaluate the non-thermal components of load. The optimization 
procedures discussed in this paper are still valid for sensor adopting this compensation 
technique, if applied to the new characteristic matrix C'. 
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